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We discuss the use of the Born and Markov approximations 
in describing the dynamics of an atom laser. In particular, 
we investigate the applicability of the quantum optical Born- 
Markov master equation for describing output coupling. We 
derive conditions based on the atomic reservoir, and atom dis- 
persion relations for when the Born-Markov approximations 
are valid and discuss parameter regimes where these approx- 
imations fail in our atom laser model. Differences between 
the standard optical laser model and the atom laser are due 
to a combination of factors, including the parameter regimes 
in which a typical atom laser would operate, the different 
reservoir state which is appropriate for atoms, and the dif- 
ferent dispersion relations between atoms and photons. We 
present results based on an exact method in the regimes in 
which the Born-Markov approximation fails. The exact solu- 
tions in some experimentally relavent parameter regimes give 
non-exponential loss of atoms from a cavity. 

42.50.Vk,03.75.Be,42.50.Ct,03.75.Fi 



I. INTRODUCTION 

An atom laser is a device which, in analogy to the op- 
tical laser, would emit a coherent beam of bosonic atoms. 
Various models of atom lasers have been presented ji]-f§ . 

Many of these schemes are based around a master 
equation description for the dynamics of the system 
|j]-f|,||. In these, the atom laser output is described by a 
Born-Markov master equation. Making the Born-Markov 
approximation involves assuming that reservoir correla- 
tions decay rapidly compared with the system damping 
time and that the reservoir does not change significantly 
with time due to the effect of the system. We discuss 
these approximations in the context of atom output cou- 
pling. 

Born-Markov master equations were initially derived 
for optical systems || . There they are used to describe a 
system (for instance an optical laser mode) coupled to a 
large, unchanging reservoir (the external modes). In op- 
tics the Born-Markov approximations allow an equation 
containing only system variables to be derived. One of 
the fundamental assumptions involved in deriving such a 
Born-Markov master equation is that the reservoir corre- 
lation functions decay rapidly. This allows the reservoir 
to be approximated as unchanging in time. While it is 
assumed the system does not affect the reservoir, the 
reservoir does affect the system. 



An equation in terms of system variables alone is also 
an important goal for describing atom lasers. However, 
atom and photon devices work in substantially different 
parameter regimes. Moreover, atoms and photons have 
different dispersion relations, which affects the behaviour 
of the reservoir correlation functions. Furthermore, typ- 
ically the reservoir for an optical cavity is taken to be 
at thermal equilibrium at a nonzero temperature. For 
atoms, a vacuum reservoir with all modes initially empty 
is often more appropriate. 

In this paper we look at the validity of the Born- 
Markov master equation for describing output coupling 
from a single mode atomic system to a reservoir. The 
reservoir is described as a continuum of free-space modes. 
Our description of the coupling will initially be quite gen- 
eral, though our later discussion will focus on particular 
schemes in which the atom becomes free of the system 
through a change of state. Such a change of state to 
an untrapped state can be achieved using either a Ra- 
man transition H or an RF transition juj-O]. We will 
also discuss the effects of gravity on output coupling. In 
the latter part of the paper we will focus on broadband 
coupling. This allows a comparison with exact results 
p5|,^6[ , however we also present finite coupling results 
which more accurately describe output coupling through 
a change of state. We will discuss exact equations which 
can be obtained in regimes where the Born-Markov ap- 
proximations fail. 

In section [ij we present our model of output coupling 
from a single mode trap. This allows exact equations of 
motion, and their solutions, to be obtained for system 
variables. In section III we model this system using the 
master equation formalism, emphasising the importance 
of the Born and Markov approximations. In section IV 
we then consider the validity of these approximations. 
In the optical case, the adequacy of the Born-Markov ap- 
proximation means that the standard model we are using 
leads to exponential decay of the number of photons in 
a cavity. The equivalent model for atoms does not lead 
to an atom number which tends to zero for sufficiently 
long time. In section VI we discuss the reasons for this 
and consider the results obtained from models which in- 
clude the effects of gravity. Having noted that the Born- 
Markov approximations may fail, we discuss in section 
VII a non-Markovian master equation and show that for 
atoms the Born-approximation is not valid for our pa- 
rameters even if we avoid making the Markov approxi- 
mation. In the regimes where these approximations fail, 
the output field becomes correlated with the system and 
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hence cannot be traced over to provide a master equation. 
Finally we discuss ways of proceeding when the Born- 
Markov master equation fails. The most straightforward 
method is that used in section [n| where exact equations 
for the whole system and reservoir are obtained. 



II. EXACT SOLUTIONS 

Dilute gas Bose Einstein Condensates (BEC) are now 
available in the laboratory. To produce a continuously 
running atom laser from a BEC requires the addition of 
a suitable pumping mechanism, and an output coupler. 
A change of atomic state through RF transition has been 
used to produce a "pulsed atom laser" |ri| , |l2|| . The out- 
put coupling from a single mode to a large reservoir is 
sometimes described by a Born-Markov master equation 
of the form 



dp(t) 
dt 



C{2ap(t)a) - a)ap{t) - p{t)a) a) 



(1) 



where C describes the strength of the coupling and a 
[a)) is the annihilation (creation) operator for the single 
mode system. Two important approximations involved 
in such a description are that the evolution is Markovian 
and that it depends only on the system operators. The 
Markovian property means that the rate of change of the 
state depends only on the state at that time. There is 
no explicit dependence on the state at previous times. 
The equation is a function of system variables only, due 
to tracing over the reservoir. This is appropriate if the 
reservoir remains uncorrelated with the system. In fact, 
it is approximated in deriving Eq. (|lj) that the reservoir 
does not change with time. 

We wish to investigate these two approximations, and 
their validity for describing atomic output coupling. In 
a full atom laser model it is essential that a pumping 
term is also included. However, in this paper we consider 
a single mode coupled only through output coupling to 
the outside world. Experimentally this corresponds to a 
leaky BEC. We focus on the output coupling term which 
would be present in a full master equation model along- 
side other terms. 

We begin by considering a generic output coupling 
mechanism which we have previously analysed in the con- 
text of Heisenberg equations of motion [jl5 16 1. Here we 
investigate how such a model can be described by a den- 
sity operator equation - in particular a master equation. 
When the output field is constrained by a waveguide such 
as in a hollow optical fibre model j^], the output field is 
effectively one dimensional. We consider a single mode 
system (the lasing mode with creation operators a)) cou- 
pled to a one-dimensional continuum of free space modes. 
The free space modes are labelled by their momentum, fik 
(creation operators b\). In a general three dimensional 
model bk would also be labelled by indicies describing 
the transverse modes. Here, however, we suppress these 



indicies and assume that the output is coupled into a 
waveguide with the atoms remaining in a single trans- 
verse mode. The Hamiltonian describing the system is 
given by 



H = Hs + Hx + H ST , 
Hs = Tiuoq a) a, 

/>oo 

H x = I dk hu> k b\b k , 



-it} 



dk («(£;) bko) — n*{k) b\a). 



(2) 
(3) 

(4) 
(5) 



The function n(k) describes the shape of the coupling 
in k space, and is left general here. By appropriately 
choosing n(k) we may simulate a wide range of practical 
output coupling mechanisms, including position depen- 
dant coupling and the effect of momentum kicks. The 
latter may result from laser photons when an atom un- 
dergoes a change of state. Choosing n{k) as constant 
over a broad region corresponds to broadband coupling. 
This model is used for optical input-output theory JlO| , 
and in proposed atom laser theories which result in Born- 
Markov master equations. For the atom case this model 
ignores potentially important effects such as gravity and 
atom-atom interactions. For now, we neglect these effects 
in order to investigate the validity of the Born-Markov 
approximations. We use this model to understand the 
differences between optical and atomic output coupling. 
We will extend the model in section (v|. 

We use the Hamiltonian presented above to write 
Heisenberg equations of motion for the operators a, bk. 
We can also obtain equations for combinations of these 
operators which may be of more interest, such as the 
number of atoms in the system, a^a. These equations can 
be difficult to solve in general. However, since they in- 
clude the output and system explicitly, they describe the 
dynamics of the model exactly |L5|,|l6|. Exact solutions 
can be compared with Born-Markov master equations. 

To facilitate this comparison we present exact equa- 
tions for the expectation values of the system operators, 
a' and a' a. In these we assume (as we will do in the mas- 
ter equation descriptions to follow) that the reservoir is 
initially empty, (6^(0)6^(0)) = 0. We do not place any 
further restrictions on (b' k (t)bk(t)) . This is the first fun- 
damental difference between the atom and photon case. 
An empty reservoir for photons is inapropriate at finite 
temperatures |]l7| . In experiments where a Bose-Einstein 
condensate is allowed to leak out of a trap, however, the 
most appropriate initial state for the outside atom modes 
is a vacuum. Similarly, for an atom laser in a hollow fi- 
bre, which we discussed in R], the initial output modes 
would be empty. 

We obtain 



dt 



(a\t)) = iu (a\t)) 



dr r(T)(aHt-r))e^ T , 



(6) 



2 



dt 



[a\t)a(t)) 



dr /'(r) {a){t)a{t-T))e- tuJ0T 



o 

h.c. 



These equations have solutions given by [[15 1G 



(o + (t)) = (a t (0))e^ ot /T 1 



a + £{/'(*)*}(*) 



(7) 



(*), (8) 



(a f (t)a(i)) = (a f (0)a(0)) x 

zr 1 



1 



s + £{/'(<)} (*) 



(*) 



(9) 



where C and £ 1 denote Laplace and inverse Laplace 
transforms respectively. f'(t) is the function defined by 



f\t) = / dfc| K (/fc)i 2 e - j( ^- wo) *. 



(10) 



The function f(t) = f'(t)e~ luJot is the "reservoir correla- 
tion function" in the master equation picture which we 
describe in the next section. Here ujk = hk 2 /(2m) for 
atoms, in contrast to u>k = CL,\k\ for photons. Here m is 
the mass of the atoms, and cl is the speed of light. 

For a previously described output coupling from a con- 
densate through state change jig] , n{k) is a Gaussian of 
width Ofc, 



K (fc)= l F^(2^)- 1 /4 e xp (-e/{Aol)). 



ill) 



Here, the strength of the coupling is given by the cou- 
pling constant, T. For this form of coupling f'(t) may be 
evaluated as: 



f(r) = 



T 



vT + iar 
where we have defined 

a = Tiaf,/m. 



(12) 



(13) 



For the broadband limit of Eq. ( |13[ ) (discussed in sec- 
tion |^) we may use methods similar to those given in 
reference |T7J] to obtain an analytical form for the inverse 
Laplace transform, 



l 



(t) 



s + £{f'(t)}(s) 
x {W(a,b,c)+W(b,a,c) + W(c,b,a)), 



W(a,b,c) 



aV * (l + erf[aVt]) 



(a — b) (a 



(14) 
(15) 



The variables, a,b and c are the three solutions of the 
equation s 3 + iuj s + Tcy/i = 0. Note that Eq. (||) gives 
that (a(t)) will always remain zero if (a(0)) = 0. For 
the case of damping of a BEC which we consider here, 



the initial state corresponds to a BEC in an atom trap. 
According to spontaneous symmetry breaking arguments 
BECs are in coherent states with a definite global phase 
so that (a(0)) ^ 0. This is a useful assumption. 
Nevertheless, even if (a) = 0, the form of the equation 
for (a(t)} must be as given. 

The equations of motion given by Eqs. (^,0) and the 
corresponding solutions, Eqs. are exact for the sys- 

tem under consideration. In specific cases it is very diffi- 
cult to solve for these inverse Laplace transforms. More- 
over, the Heisenberg equations for system operators de- 
pend on the external operators, 6^(0) in general. We next 
investigate equations of motion based on the Born and 
Markov approximations. These are compared with the 
exact solutions given above. 



III. THE BORN-MARKOV MASTER EQUATION 

Derivations of the Born-Markov master equation for 
a general system reservoir interaction are given in ref- 
erences ]l7|-pj[. Here we present a derivation for the 
specific model of Eq. (§) to Eq.(g). We assume that 
the atom reservoir is initially in a vacuum state - that is, 
there are initially no atoms outside the system. This as- 
sumption was also made in the exact solutions presented 
in Sec. |TJ. We make the Born approximation. It in- 
volves ignoring correlations which may arise between the 
system and reservoir and ignoring any time evolution of 
the reservoir density operator. We use the interaction 
Hamiltonian in Eq. (^) . This leads to the non-Markovian 
master equation: 



dp 



dr {a) ap{t — r) — ap(t — r)a) } x 



f'(r) +h.c, 



(16) 



where, p is the density operator in the interaction picture. 
The function, /'(r), defined as /'(r) = /(r) e iu ° T is the 
same as defined in Eq. ( |l0| ) and Eq. (|l2|). Here /(r) is 
the reservoir correlation function. 

Eq. (|l6|) is non-Markovian. The second major ap- 
proximation required to produce a Born-Markov master 
equation is the Markov approximation. The Markov ap- 
proximation is made on the assumption that the reservoir 
correlation function, / (t) goes to zero rapidly compared 
with the time scale on which p(t) changes. Making the 
Markov approximation thus involves replacing the terms 
p(t — t) in Eq. ([TJ]) with p(t). In the optics case, this 
approximation is also usually accompanied with the ex- 
tension of the upper limit of the integral from t to infinity. 
Making these approximations gives 

p = { (c + c*)ap(t)a) — c a)ap(t) — c*p(t)a) a] , (17) 
where 
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r t— *oo 



oo 



c= J dr I Inik^e-^-^dk 



OO 



/'(r)dr 



rV27T GXp [— [jJo/a] 



1 + Erf 



(18) 



The upper integration limit t has been extended to 
oo, as is done in the optical case, without affecting our 
subsequent conclusion regarding the Born and Markov 
approximations. This produces an equation of the form 
which has been used previously to describe atom lasers. 
We further note that we could redefine c to be real with- 
out loss of generality by incorporating the imaginary part 
of c in with the free part of the system Hamiltonian. 
This reduces the form of the loss term to the familiar 
C{2apa) — a)ap — pa^a), with C = TZe[c] the (real) cou- 
pling strength. 

The value of the constant c depends on the form of 
f'(r). In the following, we consider f'(t) to be either 
that resulting from the Gaussian coupling presented in 
this section, or the broadband limit of this function which 
we discuss in the next section. 



IV. TIMESCALE CONDITIONS 

The Born-Markov master equation, Eq. (|l7|), is the 
form used recently in various discussions of atom laser 
dynamics This Born-Markov master equation is 

used ubiquitously in quantum optics. The validity of the 
Markov approximation depends on the interplay between 
the reservoir correlation timescale, the system timescale, 
and the timescale on which the system decays. 

The condition for the validity of the Born-Markov ap- 
proximations is given in standard optics texts by the 
timescale separation condition [p~7|-|iT 



t R << At << t D 



(19) 



where tu is the reservoir correlation time, and to is the 
cavity decay time. At defines a coarse grained timescale 
on which the equations of motion are valid. Generally tu 
is defined as the "timescale on which the reservoir cor- 
relations are non-zero" . to is the decay timescale of the 
system which can be obtained by solving the equations 
for system and reservoir self consistently. In the Markov 
approximation this timescale goes as l/(c + c*) for c sim- 
ilar to that defined in Eq. ( |l8| ) except with the optics 
dispersion relation. 

Both tn and to depend on the function /'(t), and thus 
in turn on the form of LUk as a function of k. For atoms 
u>k = hk 2 /(2m), whereas for photons LUk = ci,\k\, where 
cl is the speed of light. It also depends on factors such 
as the nature of the reservoir and the parameter regime 
in which atom lasers operate. 

A second timescale condition for the Born-Markov 
approximation which is discussed in some treatements 



Pq,f9[ °f the optical Born-Markov approximation is that 
the system timescale, defined as t s = 1/wq must satisfy 



t s « At. 



(20) 



This is equivalent to requiring loq to be very large. A 
large ujq condition is required in optics for a number of 
reasons. First, the initial coupling Hamiltonian, of the 
form (ab k +a^bk), is in the rotating wave approximation 
and ignores terms of the form abk and a)b\. This rotat- 
ing wave approximation in optics can only be made in 
the case of large u>q. For the atom coupling, however, 
the correct form of the coupling does not include such 
(non atom number conserving) terms, even for small ujq 



The terms which are eliminated in the optics case [|18|]19 
through the assumption of large ujq are already zero for 
our model, due to the assumption of an atom vacuum 
reservoir. Thus, one may be led from these treatments 
to suppose that the Born-Markov approximation is made 
independently of ujq for an atom-vacuum reservoir. This 
however is not the case as we will discuss later. 

In optics, these timescale conditions are usually satis- 
fied. For a coupling based on a mirror it is standard to 
assume that the coupling is broadband [9J. That is, we 
assume n{k) is flat in k-space. In this case the reservoir 
correlation function, /'(t), given by Eq. ( |Io|) becomes 



-' luJkT dk 
- lCLkT dk 



f(r) « Mfco)|V^ / , 

= 2Kfc )| 2 / 
«2Kfc )| 2 *(t). 



(21) 



In the final equation the Dirac delta function, 5(r), is 
obtained by extending the frequency integral into physi- 
cally unrealistic negative frequencies. This is a standard 
technique in optics |J where u>q is typically large. 

Typically, for a laser, ojq « 10 15 s _1 is large and the 
assumption of a Dirac delta function decay is very good. 
Atom traps work in rather different parameter regimes 
with ujq » 10 s . If we avoid using negative frequencies, 
with the assumption of the empty reservoir considered 
here, we obtain the sum of a Dirac delta function, and 
an imaginary part corresponding to the Cauchy principal 
value of the integral in Eq. (|2l]) . 



f(r) 



2tt 



\ K (k )\ 2 e^ T ( S(r) 



CL 



I 
TXT 



(22) 



However, for an optical reservoir this estimate of corre- 
lation function decay based on our reservoir model is not 
strictly valid. This is because, while we have considered 
here the photonic dispersion relation, a vacuum is un- 
realistic for an optical reservoir at finite temperatures. 
More appropriate is a thermal reservoir, which leads to 
a decay time of order H/k B T ps 10~ 13 s 0. 

These reservoir correlation times must be short com- 
pared with the decay time, The system timescale, 
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X/ujq must also be short compared with to for a stan- 
dard optical reservoir, to is the timescale of exponential 
decay, given by e~( c+c )' in the Born-Markov limit. From 
an equivalent derivation of the optical master equation to 
that given in section H, the decay timescale is given by 
to = l/(c + c*) oc (crfe/r). That is to is inversely pro- 
portional to the strength of the coupling, given by T/at 
in the broadband limit. We will see later that for the 
atom coupling, the different dispersion relation makes to 
depend on loq and other parameters, such as the mass of 
the atom. For optical systems this decay time is typically 
of the order 10~ 8 s. Thus, in typical optical systems, the 
Born-Markov approximation holds for a number of rea- 
sons. The condition tn << to holds, because in the 
large uq limit the reservoir correlations decay as a Dirac 
delta function, to does not depend on luq and is typi- 
cally many orders of magnitude larger than the reservoir 
decay times. Similarly, the system timescale t s for real- 
istic optical cavities is very much shorter than the decay 
timescale, to- 
la contrast, the large lu limit is not necessarily valid 
for realistic atom traps. Moreover, even in the limit of 
infinitely large luq, the atom correlation function /'(r) 
does not tend towards a Dirac delta function. This is 
due to the atomic dispersion relations, which lead to a 
dependence of to on parameters other than the coupling 
strength. Furthermore, the assumption of an initially 
empty reservoir is realistic for atoms. For atoms, the 
broadband limit of the function f'(t) is given by 



/'(r)«KMI 



= 2\K(k )\ 



2 iu/ T 



2 iuiQT 



dk 



-ihk 2 /(2m)T ^ 



KM 2 



Tit 



(23) 



This is the broadband limit of the general re serv oir cor- 
relation function, Eq. (12) given in Sec. III. Both 
broadband and Gaussian coupling give forms for /'(r) 
which fall off as inverse y/r. The broadband limit of Eq. 
( |l2| ) is obtained as both and T tend to infinity, with 
r/cfc = const. = v / 27t|k(/co)| 2 - &k and V are both defined 
in Sec. H with giving the width of the Gaussian cou- 
pling, k[P). We note that the broadband limit of Eq. 
( |l2| ) is not correctly obtained by taking au — ► oo while 
keeping Y constant. If — > oo then the coupling func- 
tion n{k) and the constant c in the master equation, Eq. 
(|l7|), go to zero everywhere due to the normalisation of 
our coupling, Eq. (|ll|). 

We may now highlight three significant differences 
between optical and atomic output coupling. Firstly, 
Eq. (^) shows that the atomic reservoir correlation func- 
tion decays as \j\fr. This is different from the optical 
case even for finite ujq, and will lead to different behaviour 
of the exact equations. Secondly, the system timescale, 
is significantly larger for atom traps than for opti- 
cal cavities. Thirdly, the timescale on which the system 



decays is given by to ~ 1/ (c+c*) where c is related to the 
integral of the correlation function, /', given in Eq. (|lg|). 

The optical dispersion relation causes the system decay 
time to depend only on the coupling strength, T/a^, and 
the speed of light, c^. For the atom dispersion relation, 
to also depends on the system frequency, loq and the 
mass of the atom, m, and is given by 



to 



ugh /CTfe\ 
rmr \ T J 



(24) 



These timescale considerations based on the differ- 
ing dispersion relations, reservoir model, and parame- 
ter regimes for atoms compared with optics suggest that 
the Born-Markov approximation may not be generally 
valid in describing output from practical atom lasers. In 
fact, the optical Born-Markov theory is not universally 
valid for optics, either. In a photonic band gap material 
the dispersion relation for the photons and the radia- 
tion reservoir may be modified. Bay et al. pl| ] have dis- 
cussed fluorescence into a radiation continuum in which 
a band gap with dispersion relation near the band edge, 
LOk = uj e + A(k — fco) 2 j similar to the atomic dispersion 
relation, is present. They find behaviour which cannot 
be described using the Born-Markov theory. 

In the next section, we discuss equations of motion ob- 
tained from the Born-Markov master equation for (a^(i)) 
and (a* (t)a(t)} . We will compare these with those ob- 
tained from the exact equations, Eq. (||) and Eq. (0) 
and hence examine the regimes of validity of the Born- 
Markov approximation. 



V. THE VALIDITY OF THE BORN-MARKOV 
APPROXIMATION. 

We consider first the Born-Markov master equation, 
Eq. (|l7j). Using the relation, 



dt 



Tr 



„dp 
dt 



(25) 



where o is any system operator in the interaction picture, 
we obtain the following equations of motion for (aX) and 
tnla) from the Markovian master equation, Eq. (|17| 

w/„t/V\\ 

(26) 
(27) 



from the Markovian master equation, Eq. (|17| ) 

= [iujq-c } (a T )(£), 

d{a^(t)a(t)) f 
J t = -{c + c )(a' [t)a{t)). 



We compare the exact equations, Eq. (||) and Eq. (0) 



with Eq. and Eq. 



respectively. The equa- 



tions derived from the Born-Markov master equations 
are equivalent to the exact equations under the approx- 
imation that the term (a f (t - r))e 4 "° T = (at(t)) and 
(at(i)a(t - r))e-^° T = (a*(t)a(t)). That is, if we ig- 
nore the effect of the interaction on the system evolution. 
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Alternatively, the exact and Born-Markov equations will 
agree if /'(r) can be approximated by a Dirac delta func- 
tion. For atoms, however, there is no limit in which /'(r) 
tends to a Dirac delta function. 

When /'(t) is not given by a Dirac delta function, the 
approximation to replace (<v(t — r))e iu>aT by (aJ(t)) in 
Eq. (^) will still be valid in some parameter regimes. 
In particular, if we consider the exact equation, Eq. (]?]), 
and the solutions obtained from the Heisenberg equations 
of motion, we can see that the exact equation can be 
rewritten as 

l<ot(t)a(t)) = -<ot(t)o(t))x f f'(r) 



dt 



Jo 

X i — ; ^ ar + h.c. 



s+C{f'(t)}(s) 



}(*) 



(28) 



From this form of the exact equation, the Born-Markov 
equation is obtained by the assumption that /'(r) de- 
cays rapidly on the timescale on which the other (inverse 
Laplace transform) terms in the integral change with r. 
For parameters in which the Born-Markov approxima- 
tion is valid, we know that this ratio, given exactly from 
Eq. (|l4|), is approximately exponential with a timescale 
of order to- This fact can be motivated by considering 
Eq (j^). This equation shows that the number of atoms 
in the cavity as a function of time is given by the square 
of the absolute value of the inverse Laplace transform 
term, identical to that found in Eq. ( p8| ) above. In the 
Born-Markov regime we know the number of atoms in 
the cavity must decay approximately exponentially on 
the timescale to- Thus, for the Born-Markov approxima- 
tion to hold we require that the the timescale on which 
f'{t) decays must be small compared to to- 

For the non-broadband case, we can define a timescale 
on which f'(t) decays by the width at half maximum of 
the absolute value of the reservoir correlation function 

l/MI 



tR 



15m 



(29) 



Here m is the mass of the atoms and is the 
width of the Gaussian lasing mode in momentum space. 
This timescale must be small compared with the decay 
timescale, to discussed earlier. This condition, by itself, 
would suggest that as the coupling becomes increasingly 
broadband the Born-Markov approximations become in- 
creasingly good. However, this is not the case. Although 
the function /'(r) becomes infinite at zero in the broad- 
band limit and therefore must have a zero half width, 
tR, there are significant contributions to the integral in 
Eq. (||^) away from r = 0. Instead of the half width 
of the reservoir correlation function, tR, we are more in- 
terested in the half width of the integral of /'(r). This 
timescale is defined in terms of the reservoir correlation 
function, /(r), and the system frequency, u>q. 



We define t rs such that 

rt ra 



f'(r)dr = 1/2 / f{T)dT. 
Jo 



(30) 



For broadband coupling we find that t rs = l/u>o is 
equivalent to the system timescale, t s , defined earlier. 
The atomic dispersion relations make to also depend on 
ujq, Eq. fl24| ) which means that the condition t rs « to, 
can be written as 



3/2 



mir VI/ 



» 1. 



(31) 



For broadband coupling, this timescale condition de- 
termines the parameter regimes in which the Born- 

3/2 

Markov approximation is valid. The dependence on u> Q 
and on the mass in this condition comes from the de- 
pendence of to on m and luq. In the equivalent model 
with optical dispersion relations, to only depends on the 
strength of the coupling, given by F/ufc. 

We now compare the results of the exact and Born- 
Markov equation. We initially consider Gaussian cou- 
pling with similar parameters to those discussed in [fl6|| . 
For atoms, mw5x 10~ 26 kg. Atom traps in which BEC 
has been achieved have frequencies of uiq w 2ir x 123s -1 
p3[ . Values for the coupling strength, L depend on the 
method used. For Raman coupling, for instance, L de- 
pends on the intensity of the lasers jl6| , so a range of val- 
ues down to zero can be achieved. The value T ss 10 6 s~ 2 
can be achieved with laser intensities similar to those 
discussed in Jl6|] . The width of the Gaussian we assume 
to be determined by m 10 6 m _1 corresponding to a 
ground state wavefunction of size ~ 2/im. 

In Fig |j]a the solution for the expectation value of the 
number of atoms in the system at time t is plotted for 
the parameters quoted above. The exact solution is given 
by Eq. (||). The solution derived from the Born-Markov 
master equation is also shown. This is given by 



a T (i)a(i)) = (a T (0)a(0)> 



-(c+c*)t 



(32) 



Figure [l] demonstrates that the results for the number of 
atoms using the Born-Markov approximations disagree 
with the exact solutions in the case of Gaussian out- 
put coupling. For these parameters, tR w 2.0 x 10 _3 s, 
t s m 1.4 x 10" 3 s and t D « 5.0 x 10" 4 s, so that both the 
inequalities tR < to and t s < to fail. 
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FIG. 1. A comparison between (a^(t)a(t)) found us- 
ing the Born-Markov master equation (dashed line) and 
the exact solution (solid line). Atom numbers are nor- 
malized so that {a" 1 (0)a(0)} = 1. Parameters are (a) 
r = 1 x 10 6 s~ 2 , <Tfc = 10 6 m~ 1 for Gaussian coupling and 



(b) KM] 2 = 1.0/ 

parameters are m - 



27rms 
5 x 10" 



for broadband coupling. Other 
26 kg, luo = 2tt x 123 s" 1 . 



We have presented results in terms of the Gaussian 
coupling discussed in Sec. (III). This reflects a possible 
realistic output coupling method. However from the in- 
equality t s < tjj we see that even in the broadband limit, 
atoms do not couple out of the system in a Born-Markov 
manner for parameters which correspond to the ones we 
discuss above. Figure [Jb compares the exact and Born- 
Markov solutions to our model in the broadband limit. 
In this limit the timescale, as defined above tends to 
zero, thus the first inequality is satisfied. However, the 
inequality, Eq. (|3l]), fails. This figure demonstrates that 
even for infinitely broad atom coupling the Born-Markov 
approximation may fail. The parameters chosen are the 
same as those for the Gaussian coupling, with a strength 
of coupling in the exact solutions given by 



/ 27T (Tfc 



0.4m s~ 



(33) 



In this case the Born-Markov solution again gives expo- 
nential decay. However, the decay constant is now given 
as the broadband limit of Eq. (|l|) . 

Figure [l] demonstrates that for our model of output 
coupling from an atom laser, results for the number of 



atoms using the Born-Markov approximations disagree 
qualitatively with the exact solutions. One of the effects 
of not being able to ignore the back-action from the reser- 
voir is that for the model we are considering the number 
of atoms in the cavity does not tend to zero for long 
times (Figure [j]). We discuss reasons for this behaviour 
in the next section. If effects such as gravity and re- 
pulsive interactions are included the cavity number will 
tend to zero. However, even with these effects included 
the loss is not exponential and the conclusion that the 
Born-Markov approximation fails to describe the output 
coupling remains. 

Here, we have demonstrated the failure of the Born- 
Markov approximation by the use of the particular sys- 
tem variable, (a^a). However, the problems with using 
the Born-Markov approximations are not confined to this 
particular example. For instance, if the output from a 
BEC is described using a Born-Markov master equation, 
the resulting long time energy spectrum is Lorentzian. 
However, if we avoid making the Born-Markov approx- 
imations for atoms the exact spectrum may be non- 
Lorentzian for some values of coupling strength, T, and 
frequency, uiq Jig]. 



VI. EFFECT OF GRAVITY ON THE MODEL 

In this section we present a quasi-single particle model 
which allows us to consider the effects of gravity on our 
output coupling. Such a model cannot be extended to 
show interesting effects, such as noise suppression due 
to gain saturation if pumping is added to the model as 
it does not give information about the general quantum 
statistics of the output atom field. However, it shows 
that the inclusion of gravity causes the atom number to 
asymptote to zero. It does this in a non-exponential way, 
however, and therefore cannot be modelled by a damping 
term based on the Born-Markov approximation. More- 
over, we demonstrate that the short time behaviour pre- 
dicted by the models with gravity agree with the exact 
solutions we have presented earlier which ignore the ef- 
fects of gravity. 

The previous section demonstrated qualitative differ- 
ences between the exact solution of our model and the 
solution which uses the Born-Markov approximation. In 
fact, the exact solution of the model has a stable, nondis- 
persing state which means that not all of the atoms leave 
the cavity, whereas the approximate solution shows an 
exponential decay of atoms from the cavity. The pres- 
ence of the stable state is due to a coherence between the 
atoms in the cavity and the output modes. The Born- 
Markov approximation ignores any coherence between 
the cavity mode and the output modes, and therefore 
cannot describe any model which produces such features. 
We now show that such features would be destroyed by 
gravity. 

For coherent dynamics without atom-atom interac- 
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tions, the multiparticle evolution is identical to the evo- 
lution of a single particle 0]. The gravitational term 
makes it impossible to derive analytical results as in sec- 
tion [n], so we solved the corresponding time-dependent 
Schrodinger equation numerically. This was done in the 
position basis rather than the momentum basis for con- 
venience. The Hamiltonian for our system with the in- 
clusion of a gravitational potential is 



H = H s + H r + H sr + H g , 



H s = huj \ip a )(ip a \, 

p2 



2M' 



(34) 
(35) 

(36) 



H sr = J dx (g*(x)\ip a )(x\ + g(x)\x)(il> a \), (37) 
H g =Mgxsm(6), (38) 

where the coupling function, g(x) is related to re(fc) by a 
Fourier transform, 



g(x) = dk 



ihn*(k) e l 



/2tt 



(39) 



P is the momentum operator. With the term H g = 
0, this model is equivalent to a single particle version 
of our earlier many particle description, Eq. (|^), and 
produces a time dependance for the probability of finding 
an atom in the cavity which is identical to (a'(t)a(t)), 
found previously. These show a long time steady state in 
the number of atoms in the cavity which does not tend 
towards zero. With the inclusion of gravity {H g ^ 0), 
however, the number of atoms decays to zero in a non- 
exponential manner. This behaviour is shown in Figure 
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FIG. 2. The effect of introducing gravity into the model on 
the number of atoms in the cavity (a}(i)a(t)). Atom num- 
bers are normalized so that (a t (0)a(0)} = 1. The solid line 
shows the figure with gravity turned off and agrees with the 
results presented in Fig. |l|a. The dashed line includes gravity, 
g sin(6>), g = 9.8m s _1 ,0 = tt/20. 



In figure g we can see collapses and revivals in the num- 
ber of atoms in the cavity. This is due to the evolving 
phase relationship between the atoms in the cavity and 
the atoms which have been coupled into the output field 
modes. There is no version of the Born-Markov approxi- 
mation that can describe behaviour such as this, as such 
an approximation requires that there be no entanglement 
between the lasing mode and the output modes. 

The presence of gravity changes the long time be- 
haviour of our model so that the number of atoms asymp- 
totes to zero, while the short time behaviour remains the 
same. Other effects may also lead to the long time decay 
of atom number. Repulsive interactions, for instance, 
would be expected to destroy the thin stable structure 
in position space which leads to the long time non-zero 
population of the cavity mode. The effect of repulsive 
interactions can be included in this model by including 
in Eq. (p4j) a nonlinear Hamiltonian term given by 



Hi = NUoW(x)\ 2 



(40) 



Including such interactions into our model produces a 
Gross-Pitaevskii type equation ^0|, where N is the 
number of atoms in the system and Uq an interaction 
strength. We find that introducing such an interac- 
tion term does cause the atom number to go below the 
nonzero steady state predicted in the interaction free 
model. 



VII. A NON MARKOVIAN MASTER EQUATION 

In the previous two sections we have shown that the 
standard master equation does not correctly describe the 
dynamics of our model for particular parameters. How- 
ever a master equation, in terms of only the system vari- 
ables, would be an important tool for describing an atom 
laser. We now consider whether a non-Markovian master 
equation can give a correct description of the atom laser. 
To do this, we continue to make the Born approximation, 
but do not make a Markov approximation. 

The master equation with the Born approximation 
only is given in Eq. (|l6|). Again, we check the valid- 
ity of the Born approximation by comparing the results 
obtained from this Born master equation with the exact 
solutions. We begin by considering the resulting equation 



for the expectation value 



^M= iwo < t)(t)- f dt'{a\t'))x 
at J Q 

{Kik^e^-^dk, 



(41) 



Eq. ( fill ) is the same as that obtained through the full 
system plus reservoir equations given in Eq. (^). This 
can be seen by making the transformation r = t — t' in 
Eq. ( f4l| ) to obtain the alternative form, Eq. ([j|). Despite 
this success, the density operator equation with the Born 
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approximation, Eq. ( Jig ) is not correct. In particular, 
the Born approximation does not give correct values for 
higher order expectation values, such as (aJa). 

The equation derived from the non-Markovian master 
equation, Eq. (|l6|) for the number operator expectation 
value is 



d{at(t)a(t)) 



dt 



= — j dr{a i a)(t — r) x 



|K(fc)| 2 e i( "°- LJfe)T dfc + h.c, 



(42) 



with solution 

(a^t)a(t)} = (a f (0)a(0))x 

c- 1 



£{/'(*)} (s) + £{/'(*)*} (s 



(*)• (43) 



These do not agree with the corresponding exact equa- 
tions given in Eq. (0) and Eq. (Q), as is shown in Fig. 
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FIG. 3. A comparison between (aJ(t)a(t)) found using the 
Born (non-Markovian) master equation (dashed line) and 
the exact solution (solid line) respectively in the broad- 
band regime. Atom numbers are normalized so that 
(a^(0)a(0)) = 1. Parameters are the same as in Fig. |l]b. 

This figure compares the exact results for (a^(t)a(t)) 
with the solution in the Born approximation only, Eq. 
( fi"3| ) . The parameter values chosen are the same as those 
discussed in section ^ and the comparison is made for 
broadband coupling. The reason for the failure of the 
non-Markovian master equation is that where the correct 
equations lead to the term (a? (t)e~ luJQT a(t — r)), the non- 
Markovian master equation with the Born approximation 
leads to {o?{t — r)a(t — r)). These two terms are equal 
if we consider only the free evolution of the system, and 
ignore the interaction term. That is, to zeroth order in 
Kfc)| 2 , (ot(t - r)a(t - r)) m (at(t)a(t - r))(i) e -^» t . 
However, there are interaction terms between the system 
and reservoir. Due to the interaction, terms obtained 
using the Born approximation are not equal to the exact 
results. To reproduce equations for the correct system 



dynamics we cannot ignore the correlations which arise 
between the system and reservoir. 

In optics the Born-Markov approximation is a useful 
approximation in almost all practical parameter regimes. 
In fact, in optics the value of luq is sufficiently large that, 
if one assumes a linear dispersion relation the reservoir 
correlation can be well approximatied by a Dirac delta 
function. This is not the case for atoms. 

The non-Markovian master equation, Eq. (]!(]), gives 
correct results for single operator expectation values. 
Nevertheless, the Born approximation neglects higher or- 
ders in the interaction term, |«;(fc)| 2 , thus ignoring reser- 
voir evolution. This leads to incorrect results for higher 
order expectation values, such as (a^(t)a(t)). 



VIII. CONCLUSION 

We have discussed the use of master equations and 
other density operator equations for atom lasers. We 
find that the Born-Markov master equation, as is com- 
monly used in optics, is not valid in parameter regimes 
in which atom lasers may work. The output modes cor- 
relation function for photons is well approximated by a 
delta function in the broadband limit, and for typical op- 
tical parameters. In comparison, this function falls off as 
1/Vt for atoms because of the atom dispersion relations. 

The Born-Markov approximations must be made self 
consistently. In regimes where the Born-Markov approx- 
imations fail, the system can be solved in the Heisen- 
berg picture, treating the output modes fully OJlq ]. For 
broadband coupling, the parameter regimes in which the 
Born-Markov approximation is valid is given in Eq. (|3~l|). 
Another possibility for describing such systems involves 
the use of quantum trajectories. These have been found 
to be useful for other systems in which no Born-Markov 
master equation exists [^lj . 

The breakdown of the Born and Markov approxima- 
tions means that new theoretical methods will be re- 
quired to understand the atom laser in some regimes. 
However, it also opens up the possibility of finding new 
properties of atom lasers significantly different from those 
found in the optical laser. 
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